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We compute the renormalized expectation value of the square of a massless, conformally coupled, 
quantum scalar field on the brane of a higher-dimensional black hole. Working in the AADD brane- 
world scenario, the extra dimensions are flat and we assume that the compactification radius is large 
compared with the size of the black hole. The four-dimensional on-brane metric corresponds to a 
slice through a higher-dimensional Schwarzschild-Tangherlini black hole geometry and depends on 
the number of bulk space-time dimensions. The quantum scalar field is in a thermal state at the 
Hawking temperature. An exact, closed-form expression is derived for the renormalized expectation 
value of the square of the quantum scalar held on the event horizon of the black hole. Outside the 
event horizon, this renormalized expectation value is computed numerically. The answer depends on 
the number of bulk space-time dimensions, with a magnitude which increases rapidly as the number 
of bulk space-time dimensions increases. 

PACS numbers: 04.62.-|-v, 04.70.Dy, 04.50.Gh 

I. INTRODUCTION 

Higher-dimensional brane-world models [in have the 
consequence that the energy scale of quantum gravity 
may be many orders of magnitude smaller than the usual 
Planck scale 10^® GeV, and may be as low as the TeV- 
scale. This raises the exciting possibility of probing quan¬ 
tum gravity effects in high-energy collisions, either at 
the LHC or in cosmic rays Of the possible quan¬ 
tum gravity processes, the creation of microscopic black 
holes would be particularly spectacular 

If such a microscopic, higher-dimensional, black hole 
were created, it will be short-lived, decaying rapidly due 
to the emission of Hawking radiation. It is expected 
that for a significant proportion of the evolution of the 
black hole it can be modelled semi-classically Q . In this 
regime the geometry of the black hole is regarded as clas¬ 
sical, with quantum fields propagating on a background 
space-time metric. The details of the Hawking radiation 
emitted by the black hole as it decays are of particu¬ 
lar importance for the simulation of black hole events 
at the LHC [l3, • There is now a vast literature on 

Hawking radiation from higher-dimensional black holes, 
see for some reviews on this subject. 

The fluxes of energy and angular momentum in the 
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Hawking radiation are just two components of the ex¬ 
pectation value of the stress-energy tensor for a 

quantum field on a black hole background. The other 
components of (T^v) also contain physical information, 
and the stress-energy tensor as a whole governs the back- 
reaction of the quantum field on the space-time geometry 
via the semi-classical Einstein equations 

G^,.=8^G(V). (I.l) 

The stress-energy tensor involves products of a quan¬ 
tum field operator at the same space-time point and 
therefore, in general, requires renormalization. The black 
holes of phenomenological interest for high-energy colli¬ 
sions either have a single axis of rotation which lies in the 
brane or are non-rotating. In both these cases, the sym¬ 
metries of the brane black hole space-time mean that the 
components of (Tjjj/) which comprise the Hawking fluxes 
of energy and angular momentum on the brane do not 
require renormalization (2ll - [^ . 

Even in four space-time dimensions, computing the 
full renormalized expectation value of the stress-energy 
tensor for a quantum field on a back gro und black hole 
space-time is a complicated process In more 

than four space-time dimensions, Decanini and Folacci 
[s^ have developed a general formalism for the com¬ 
putation of (T^j/)ren for a quantum scalar field (f based 
on Hadamard renormalization. Some general properties 
of (T)ii/)ren for higher-dimensional black hole space-times 
are derived in [s^, but to date no computation of {Tfiv)ren 
on the full space-time exterior to a higher-dimensional 
black hole event horizon has been attempted. A quan- 
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turn scalar field is mathematically the simplest type of 
quantum field, and in this case the renormalized vac¬ 
uum polarization (0^)ren is considerably easier to calcu¬ 
late and shares some features with the full renormalized 
stress-energy tensor. The vacuum polarization is also of 
physical significance for spontaneous symmetry breaking 
(see, for example, for a discussion of this in the black 
hole context). 

The vacuum polarization has been computed on a 
wide variety of black hole space-times in four dimensions 
[4(il - l5lj |. For higher-dimensional black holes, as well as 
the approach of based on Hadamard renormaliza¬ 
tion, a general formalism for computing (^^)ren based on 
de Witt-Schwinger renormalization has been developed 
[^ . An exact expression for the renormalized vacuum 
polarization on the event horizon of a five-dimensional 
asymptotically flat black hole has been found in . As 
far as we are aware, the only complete computation of 
the renormalized vacuum polarization everywhere out¬ 
side the event horizon of a higher-dimensional black hole 
is the work of on a five-dimensional, asymptotically 
anti-de Sitter black hole. 

In this paper we study the vacuum polarization (<()^)ren 
of a massless, conformally coupled, quantum scalar field 
on the brane of a higher-dimensional black hole in the 
context of the AADD brane-world scenario EH- In this 
set-up, our universe is a four-dimensional brane in a D- 
dimensional bulk space-time (where D > A). The D — A 
extra dimensions are flat but compactified. The radius 
of compactification is sufficiently small to avoid contra¬ 
dictions with experimental searches for deviations from 
Newton’s Law of Gravitation, but is typically large com¬ 
pared with the Planck length. We use a very simple 
model of a black hole in this scenario, assuming that 
the brane is tensionless and its thickness negligible. We 
also assume that the size of the black hole is very small 
compared to the compactification radius of the extra di¬ 
mensions. With these assumptions, the black hole can 
be modelled as a higher-dimensional, asymptotically flat, 
solution of the vacuum Einstein equations. For simplic¬ 
ity, we restrict our attention to static, spherically sym¬ 
metric, higher-dimensional black holes described by the 
Schwarzschild-Tangherlini metric . 

In the AADD brane-world scenario, the particles and 
forces of the Standard Model are constrained to live on 
the brane, in order to avoid contradictions with precision 
particle-physics experiments. Only gravitational degrees 
of freedom can propagate in the bulk extra dimensions. 
While it is possible to have scalar fields in the gravi¬ 
tational sector of the theory and therefore in the bulk 
space-time, studies of Hawking radiation have revealed 
that the emission of scalar fields in the bulk is suppressed 
relative to the emission on the brane . With 

this in mind, in this paper we consider only the vac¬ 
uum polarization on the brane. This means that we can 
restrict our attention to the four-dimensional on-brane 
metric describing the black hole, and therefore use es¬ 
tablished methodology [H, |4l|, [5l| to compute ((^^)ren- 


The outline of this paper is as follows. In Sec. Ill Al we 
outline the metric of the black hole on the brane, and con¬ 
struct the point-split Euclidean Green’s function which 
will be used to calculate the vacuum polarization when 
the quantum scalar field is in the Hartle-Hawking state 
[^ . We use two different methods for finding (0^)ren 
outside and on the event horizon of the brane black hole. 
In Sec. Ill Bl we follow the approach of [1^, [5l|, using 
time-like point-splitting, to find an expression for the 
vacuum polarization outside the event horizon. The vac¬ 
uum polarization on the event horizon is calculated in 
Sec. Ill Cl using radial point-splitting [4l|. On the horizon 
we have a closed-form expression for ((()^)ren, but away 
from the horizon numerical computation is required. The 
results of this computation are presented in Sec. iniifor 
D = A ..., 11 and our conclusions are in Sec. El 

II. METHODOLOGY 
A. Set-up 

We start with the metric for a D-dimensional 
Schwarzschild-Tangherlini black hole: 

= —f{r) dt^ + f{r)~^dr^ + dflD -2 (2.1) 

where 

/(^) = 1- (v)''”' (2-2) 

and dQ,D -2 is the line element on the {D — 2)-sphere 
Here and throughout this paper we use Lorentzian 
metric signature -I-) and units in which SttG = 

c = h = ks = ^- The metric (EU) represents a static, 
spherically symmetric, black hole with mass 

M=l{D-2)r^-^AD-2 (2.3) 

where Ad -2 = 27r(^“^i/^/r[(D — l)/2], with F the Euler 
gamma function, is the area of the {D — 2)-sphere. 

Labelling the angular co-ordinates on the {D — 2)- 
sphere above by 0, </?, 63 ,... , 9 d- 2 , the metric of the 
brane black hole is constructed by setting 9i = -kI 2 for 
* = 3,..., D — 2, giving the on-brane metric 

ds2 =-fir) df + f{r)-^dr^ + sin^ 6 dip^] . 

(2-4) 

In the four-dimensional brane metric (1^ . the function 
/(r) still has the form (12.2L and so depends on the num¬ 
ber of bulk space-time dimensions D. The brane metric 
(EH) is therefore not a solution of the vacuum Einstein 
equations except in the case D = 4, when it reduces to 
the usual Schwarzschild form. When D > A, the metric 
(12.41) is a solution of the nonvacuum Einstein equations 
with an effective fluid source on the brane [s^ . 

In this paper we consider a massless quantum scalar 
field (j) with conformal coupling to the four-dimensional 
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space-time geometry (EH). The scalar field satisfies the 
Klein-Gordon equation 


[V^V^ -^R]$ = 0, (2.5) 


where ^ = 1/6 is the coupling constant for conformal cou¬ 
pling on the four-dimensional brane black hole space-time 
(12.41) and R is the Ricci scalar of the four-dimensional on- 
brane metric (1^ . which takes the form 


R = 


{D - 4){D - 5)r^-^ 

j.D -1 


( 2 . 6 ) 


The Ricci scalar therefore vanishes for D = 4, 5. 

To compute the renormalized vacuum polarization, we 
follow [ 23 , [ 2 ^, [5l| and take a Euclidean approach. Defin¬ 
ing Euclidean time r = it, the brane metric (12.41) becomes 


where Puji{r) and quj(.{r) are solutions of the radial equa¬ 
tion 



^2f . 

df\ 

dG,^t 

dr'^ \ r 

dr) 

dr 

'uP- £{£+!) 

R 
“h — 

G^t 

L / 


6 J 


(2.14) 


which is the homogeneous version of the radial equation 
for arising from separating the inhomogeneous Klein- 
Gordon equation (12.81) on the metric (12.71) . The function 
Pujiir) is defined as the solution of (I2.14|) which is regular 
at the event horizon r = rh, while gujiir) is the solution 
of (12.141) which is regular as r —>• 00 . In (I2.13L as usual, 
r< is the smaller of the two values r, r', while r> is 
the greater. The normalization constant C^je in (12.131) is 
determined by the normalization condition [24 


ds^ = f{r)dT^+f{r)-^dr^+r^ [dO^ + sin^ 6 dip^] . (2.7) 

The point-split Euclidean Green’s function Ge{x; x') sat¬ 
isfies the equation Hill 

m Ge{x-x') = -g-G^{x)6\x, x') (2.8) 

where the covariant derivative is taken with respect to the 
Euclidean metric (IQ) . which has determinant g. The 
unrenormalized expectation value of the vacuum polar¬ 
ization (/>^)unren IS given by the coincidence limit 


(</' 


")unren = K 


lim Ge{x; x) 


(2.9) 


We consider the quantum scalar field to be in the ther¬ 
mal Hartle-Hawking state at temperature T. The 
four-dimensional black hole (12.41) has Hawking tempera¬ 
ture 


T = 


D-3 

Aeth 


( 2 . 10 ) 


which increases linearly with D for fixed Vh, and is in¬ 
versely proportional to for fixed D. For a thermal 
state, the point-split scalar Euclidean Green’s function 
Ge(x\ is periodic in t — t' with period T~^, and takes 
the form [IJ [H, |4^ 


G^i 


Pu>e 


dquji 

dr 


- qui 


dpuii 

dr 



(2.15) 


Before we can bring the space-time points together in 
(12.91) by taking the limit x ^ x', we need to subtract di¬ 
vergent terms. For a massless, conformally coupled scalar 
field, these take the simple form 


(</')div = 


1 

87r^tT 


1 Rai3(T°'a^ 

967r^ a 


(2.16) 


where a is one-half the square of the geodesic distance 
between the points x and x' and (t“ = with R^p the 
Ricci tensor. The detailed form of the divergent terms 
(12.161) depends on the point-splitting chosen. 

For the remainder of this section, our methodology 
for computing the renormalized vacuum polarization falls 
into two parts: (a) the computation away from the event 
horizon, using temporal point-splitting following [^ . 
and (b) the computation on the event horizon, using ra¬ 
dial point-splitting following [4l| . In the rest of this paper 
we set the event horizon radius r/j = 1, so that all lengths 
having numerical values are in units of the event horizon 
radius. 


B. Outside the event horizon 


rp 

G_e(x; x') = — exp [iuj (r - r')] 

47r ^' 

n— — oo 

00 

X (2£-b l)P^ (cosy) G,^£(r;r') (2.11) 
1=0 

where uj = 2m:T, and is the usual Legendre function 
with 

cosy = cos 9 cos 6 ' + sin0 sin0' cos {ip — p'). (2.12) 

The radial Green’s function G^i takes the form 

(2.13) 


Outside the event horizon, we follow the standard 
methodology developed in [ 2 ^ to compute the renormal¬ 
ized vacuum polarization (0^)ren- We choose temporal 
point-splitting, setting r = r', 9 = 9' and p = p'. Then 
cosy = 1 (12.121) and using the fact that Pi{l) = 1, the 
Euclidean Green’s function ()2.11l) takes the form 

Ge{t,x;t',x') = 

rp CO 00 

_ ^ e*-^y^(2£+l)G^,p,,,(r)g,,,(r), (2.17) 

n— — oo ^—0 

where e = t — t' . As is well known, the sums over £ 
in (12.171) do not converge although, by definition, the 


G,,;^(r,r') = Ga,t.Pu>t{r<)q^i{r>) 
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Euclidean Green’s function Ge must be finite when the 
points are separated. This is remedied by subtracting 
from (12.1711 a suitable multiple of the Dirac delta func¬ 
tion, which vanishes when the points do not coincide, 
giving 


mode functions. We define a new function Cuitir) by 

= C^ip^e{r)que{r), ( 2 . 22 ) 

which satisfies the following differential equation [sTj l 


Ge{t,x\t' ,x') 

rri OO OO 

n— — cc> l—Q 


+ l) 


1 

'T'/f. ' 
(2.18) 


For temporal point-split ting , the divergent terms (12.1611 
take the following form for a massless, conformally 
coupled, scalar field: 


(^^)div = 


1 

dTT^/e^ 


1 

1927r2/ 


1 df 

dSTT^r dr 



1 

967r2 dr^ 

(2.19) 


These divergent terms can be written as mode sums by 
using standard identities [1^ , subtracted from p.l8[l and 
then the limit e —>■ 0 taken. Th is g ives the final renor¬ 
malized expectation value to be |25l | 


)ren — )analytic “b )numeric (2.20) 


where 


T 2 


(0 )analytic — lor 


1 


df 


12 f 1927r2/ \dr 

1 df 

dSTT^r dr ’ 


J_fl 

967r2 


-b 


(2.21a) 


T 

27T 


numeric — 

OO ( OO r 


EE 

n=l 1^=0 

rri OO 

-E 

47r 2^ 


e=o 


{21 -I-1) Cujepujeir)quje{f) - 

{2i + 1) GoePoe{r)qoeir) - 


1 UJ 

1 


Vf\ 


(2.21b) 


The mode functions Pu:i{r), quji{r) cannot be found 
in closed form for w > 0, so that, as the name sug¬ 
gests, (9l)^)numeric (|2.21bll is Computed numerically while 
we have a closed form expression for (i/)^) analytic (|2.21a|l . 

To speed up the convergence of the mode sums in 
(l23Tbll . we follow M and subtract a WKB-like approx¬ 
imation, which encodes the large w, i behaviour of the 



■ 

where 



d (^ 2 r d{V^ \ 
dr y dr J 



,(2.23) 


Xuje{r) = y a;2r-4 + 0 ^2/, (2.24a) 

V{r) = (2.24b) 


and the Ricci scalar R is given by (1^ . The WKB ex¬ 
pansion is found by inserting a fictitious parameter e in 


(IT^ : 





1 - 


1 


r'^f d ( 

7CTr 


V 


(2.25) 


expanding Ci^i{r) in inverse powers of e 

= C,0uj(.{r)+e~‘^C,iuji{r) + e~^(^2^(,{r) + e~*^C,^^i{r) + . .., 

(2.26) 

and then setting e = 1 at the end of the calculation. The 
WKB terms Ciinf(^) have the form 


2i-|-l 

E 


Qu,i{r) = E Ai^k{uJ,r)x, 


-{2i+2k-l) 

uii 


(2.27) 


where the functions Ai^k{^,r) depend on cj (and r) but 
not The first of these functions is blnjfa;, r) = 1/2 and 
the forms of ^i,fe(w, r) can be found in [^. To reduce the 
number of mode functions Puii{r), quii{r) which have to 
be found numerically, we used an expansion ()2.26p up to 
and including C,d.uji{r)\ the remaining functions A 2 ^k{'-^^r) 
and bl 3 _fe(a;,r) can be found in [^ . 

Adding and subtracting the WKB expansion from 
(l23Tbll . we arrive at the expression 
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T 

)numeric = ^ ^ ^ ^ (2£ + l) (^) Coa;^(^) Cluj^ij'') C2u)i{j'^ C3a;^(^)] 


GO I OO 


27r 


n—l k ^=0 


^(2£+l) 
f=o 
T 


Couj^i'/') + Ciw^(^) + C2u>t.{f) + Csuiei'r) — 


v7j / 


^ [^oePoeir)qoe{r) — Cooiir) — (loe(r) — ( 2 oe(r) — C30^(?')] + Ai + A 2 + A 3 , (2.28) 


^=0 


where 


Ai - ^ (2i+ l)Cio£(r) - -—jAi^i(0,r), 

r'^/ 2 


(=0 


A 2 = ^( 2 £ + l)C 2 o,(r) = 


e=o 

00 


3r^f' 

27r6 


:242,i(0,r), 


A 3 = {2i+l)CMr) = —^^ 3 , 1 ( 0 , r). (2.29) 

15r'j2 


e=o 


The mode sums in the first and third lines of (I2.28|) are 
now rapidly converging. It remains to find the sums over 
the WKB expansions in the second line of p.28|) . To 
do this, we follow [H, El, [sij and employ the Watson- 
Somerfeld identity 


Y,He) = 


i=0 


A=0 

-JR 


^ ( A - - ) dA 




/A=0 


(2.30) 


which is valid for any function J-{i) analytic in the right- 
hand half-plane. Using (12.301) we write (for w > 0) 


^(2£+l) 


e=o 


Couj^ir) - 


'V7J 


= Io{uj,r) + Jo(w,r) 
1 


24wr2 ’ 

00 

{21 + 1 ) (ku,i{r) = Ik{uJ, r) + Jk{uj, r), 

(2.31) 


e=o 


for k = 1,2,3, where the Ik{ui,r) are the first integrals 
coming from (12.301) and the Jk{oJ,r) are the second. 

The Ik{uJ,r) integrals are readily computed; the first 
two take the following form for a massless, conformally 


coupled scalar field 


lo 

h 


{^,r) = 


2 XC,ou>i — 




2A^icj^ dA — 


dA = -^, 


1 


2 Au)r‘^ 


(2.32) 


The remaining l 2 {u!,r) and l 3 (u;,r) are straightforward 
to find, for example, in Mathematica but the resulting 
expressions are sufficiently long that we do not include 
them here [^ . 

To compute the Jk integrals, we first make a change of 
variables 1 ^: 


ujr 

X = pq, ^^" 77 ’ 


(2.33) 


The integral Jq is analysed in M- After an integration 
by parts, it takes the form 


Jo = %- 


1 


Attuj 


f 24ojr^ f 


/A=0 


1 - 


A 


2 \ 2 


„27rA 


-^dA 

(1 -^g277A)2 

(2.34) 


which is 0(w ^) as a; —>■ 00 . For the remaining Jk inte¬ 
grals, using (I2.27|) . we have 


2fc-|-l 

i=i 

r g(i-g^)~ ^ 

4=0 1 + 


xJR 


dq 


for fe = 1, 2, 3 where 

y=( 2 j + 2 fc-l)/ 2 . 


(2.35) 


(2.36) 


In M , the integral over q in (12.351) was performed by re¬ 
peated integrations by parts. Here we take an alternative 
approach, first defining 


g(g) = 


1 + e'^^pq 


(2.37) 


and letting gy{q) be the Taylor series expansion of g{q) 
about g = 1 up to order y + (1/2). Then, for k = 1,2,3, 
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f/*®® 1') 

Jfc=^V(pVV) ^ Ak,jiuj,r) I [g{q) - gy{q)] (1 - q^) ^ dg + JR / gj^(g) (l - ^ dq (2.38) 

. j Ld</=o U9=o J J 


The integrand in the first integral in (12.381) is regular at q = 1 and therefore the first integral is straightforward to 
compute numerically. The integrand in the second integral in (12.3811 has a pole at q = 1. Integrating around the pole 
using the contour shown in Fig. 1 in [5l| is routine in, for example, Mathematica as gy{q) is a polynomial in q, and 
yields a finite answer in each case. 

Finally, ((^^)numeric defined by Eq. (12.281) takes the form 

J, OO oo 

(0 )numeric “ ^ ^ l ^ ^ (2^ “t“ 1) [Cuj^Pujl^T^quj^^T'j Coa;^(^) C2aj£(^) C3a;^(^)] 

n=l ^ e=o 

+ I2 + . I3 + >^0 + >^1 + <^2 + daj ’ 

rj^ 00 

+ ^ {C^o£Poe{r)qo£{r) — (oo£(r) — Ci 0 £(£') — C 20 £(£') — C30£(r)] + Ai + A 2 + A 3 . (2.39) 

^ £=0 

The results of numerically computing (12.391) will be discussed in Sec. IIIICl 


C. On the horizon 

The method discussed in the previous subsection works 
well for any r > 1 but is not well-suited to a computation 
of ((/>^)ren On the horizon r = 1 (for example, there are 
terms in the A^ (12.291) which diverge as r —>■ 1 and / —^ 
0). We therefore take a different approach to computing 
(<^^)ren On the horizon, following [ilj and using radial 
instead of temporal point-splitting. 

Separating the space-time points along the radial di¬ 
rection, we may express the unrenormalized vacuum po¬ 
larization (12.91) in the form 


The advantage of using radial point separation is now 
clear as we have a sum which involves only the zero fre¬ 
quency modes. In fact, as argued in [^ . we may only 
choose radial or angular separation for on-horizon calcu¬ 
lations as temporal separation is meaningless there. 

We renormalize (12.431) using the method of Brown and 
Ottewill , by subtracting off the singular terms in the 
Hadamard expansion of the Euclidean Green’s function, 
which we denote by ((/>^)div These have been worked 
out explicitly for radial point separation in [4l|, and for 
a conformally coupled massless scalar field (with /(r) 
given by (12.21) with rh = 1) take the form 




lim 

r—¥r' 


unren — 

00 00 


(2£ -I- 1) Cu,£Pu£{r<)qu>£{r>) 


_n=0 £=0 

where r< = min(r, r'), r> = max(r, r'), and 


F{n) = 


T/Att, 

T/ 27 T, 


n = 0, 

n > 0. 


, (2.40) 


(2.41) 


We now set = r, place the inner point, r', on the 
horizon and consider the limit r —> I. In a neighbour¬ 
hood of r' = 1, the solution p^f(r') of the radial equation 
(12.141) which is regular at the horizon possesses the series 
expansion [4lj 


PLj£{r') = 


^ (r'-l)^+o[(r'-l)^+^ 




(2.42) 


{2^)div = 


£1-3 


D-S 


167r2(r — 1) 487r^ 


-I- O [(r — 1) ln(r — 1)] . 

(2.44) 

We then have a formal expression for (^^)i.en, given by 


)ren 


lim 

r—)-l 



00 


E 

£=0 


(2£ + 1) Coiqot(r) 


D-3 D-3 

167r^(r — 1) 487r^ 


(2.45) 


We now need to calculate the sum over the modes 
qo£{r)- We begin this process by noting that the radial 
equation (12.141) . for the n = 0 modes, possesses closed- 
form solutions in terms of hypergeometric functions. To 
demonstrate this, we define new independent and depen¬ 
dent variables as follows: 


where w = 2TTnT. Setting r' = 1 therefore has the effect a; = 1 — , Z{x) = (1 — x)°‘Go£ 

that Paj£{r') vanishes for n > 0. Therefore (12.401) reduces 

to: with 


/ rin 

)unren teE {2i+l)Co£qo£{r). (2.43) 

£=0 


3(4 -D)- y/S{D-A){D-2) 
6(£)-3) 


(2.46) 


(2.47) 
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Under this transformation, (12.141) takes the form 
0 = a;(l — x)Z"{x) 


1 - 2 


Zl-3 


y/i{D - 3) 

£+1 


D-3 


-a]Z{x). (2.48) 


We may express (12.481) in the form of the hypergeometric 


differential equation [64| . 

x(l — x)Z''{x) + [c — (a + b + l)x] Z'{x) — abZ{x) = 0, 

(2.49) 

by making the following identifications 


a = — 


D-3 


b = - a, 

D-3 


c = 1. (2.50) 


It then follows that (I2.48|) has the following pair of linearly independent solutions 


Zi{x) = 2Fi{- 


D-3 


£ + 1 

— C(, — -- — a; l:x 


Z2{x) = (1-x)-^+“2Fi 


D-3 
t + D-3 
D-3 


, J. , c«y I , 


^+l 

D-3 


2i + D-2 


D-3 ’ 1-x 


(2.51) 


Transforming back to our original variables, and identifying the solutions regular at r = 1 and r = oo, we find that 
the two independent solutions to (I2.14p . poi(r) and qoe(r), are given by: 




D-3 
e+i 




qoi{r) = r 


- ^-(^+1) 


2V^r 




D-3 


D-3 


D-3 


a. 


D-3 


- a; ■ 


2£ + D-2_ 
D-3 ’ 


„3-£> 


(2.52) 


The numerical factors in front of poi{r) and qoi{r) are 
chosen so that their series solutions are in agreement 
with the expressions in (^ . but this is purely a mat¬ 
ter of convention. We could alternatively have chosen to 
include these factors into the constant Cqi (which in our 
convention is equal to unity). 

Returning now to (12.451) . while we have a closed form 
expression for qoeii^) given above (I2.52L unfortunately 
we are unable to immediately perform the mode sum as, 
to the best of our knowledge, no closed-form expression 
for this sum is known. We therefore adapt the method 
developed in [4l| to the case at hand. The procedure 
in [ilj can be applied to calculate ((^^)ren for a general 
spherically symmetric black hole space-time (I2J1) . where 
the only constraint on the function f{r) is that it has a 
single zero at the horizon in question, as is the case here. 

We begin by writing the full solution for qoe{i^) in terms 
of an approximation in the following manner 


qoe{r) = Qoe(r) + l3iPoi{r) + 77.^ (r), (2.53) 


where Qoi{r) is some approximation to qoi{r), Pi is a 
function of £ only and TZi{r) denotes the remainder of 
qnf.(r) not captured by the first two terms. The method of 
|4ll) makes use of a uniform approximation, developed us¬ 
ing extended Green-Liouville asymptotic analysis. This 
gives Qoi{r) to be 


Qoe{r) 



1/4 


Kn 


koe^Hr) 


(2.54) 


where 

fco = '\/{£ Cl) {£ -\- a) 

with 



(2.55) 


(2.56) 


for the case in hand (see [4l| for the general form), Kq is 
a modified Bessel function of the second kind and 


^(r) = 




{D - 3)- 


. In (r^D-3)/2 ^ 



(2.57) 


In (4ll | it was shown that, for a general /(r) (under 
the constraints mentioned previously) the approximation 
(12.541) encapsulates enough of the near horizon behaviour 
of qaiif) in order to capture all of the local contribution 
(i.e. the horizon divergence and local finite terms) to the 
mode sum in ()2.45l) in the limit r —^ 1. In other words, 
the sum of the remainder term TZi{r) does not contribute 
in this limit. The Pi terms in (12.531) are determined by 
the requirement that qoi(r) should vanish as r —>■ oo. 

























































In the limit r —>■ 1 we then have 


on the right hand side of (12.581) and combine them with 
the renormalization subtraction terms, before taking the 
r —>■ 1 limit, namely 


oo 


2'kT 




i=0 


1=0 


+0 [(r — 1 ) ln(r — 1 )] . 


(2.58) 


Therefore in order to calculate (^^)ren for the space-time 
(lO) . we are required to calculate the two summations 


J 


)ren — lim 
r—>-1 


^ J2(2i+l)Qo,- 


327r3 


£1-3 


487r^ 


e=o 

1 


167r2(r — 1) 


+ (2.59) 


e=o 


A closed form expression for the first term on the right hand side of (12.591) has been calculated in [41|, which here 
reduces to 


1 

9671^ 


1 2D - 5 - 2iV3hi 

[r(a )1 

\ 

[r(a)J 


+ 12-C(x,a) 



+ 12-C(x,a) 



(2.60) 


where ^{x,a) is the generalised Riemann zeta function. 

All that remains now is to calculate the /3^ summation in (12.591) . In general, this must be performed numerically 
(indeed in general Pi itself must be obtained numerically for each £), however in this case, primarily because we have 
an expression (12.521) for qoi{r) in terms of hypergeometric functions, we may calculate this sum analytically to give a 
completely closed form expression for ((/)^)ren on the brane for an arbitrary number of bulk space-time dimensions. 

Our first step in this process is to find the Pi themselves. We compare the lowest order terms in the expansions of 
the two expressions (12.521) and (12.531) for qoe{r). Expanding (j2.52l) gives: 

Qoiir) = |27 + In [(D - 3)(r + ^ + ^|; | + O [(r - 1) ln(r - 1)], 

( 2 . 61 ) 

where 7 is Euler’s constant and iji is the digamma function, while by expanding (12.531) we obtain: 

<loe{r) = |27 + In + In [(^ + a){i -f a)]| -b + 0[{r- 1 ) ln(r - 1)]. (2.62) 

Comparing the two expansions (12.611 we arrive at an expression for Pf. 

= \ jin [{(■ + a)[£ + a)] - 2 ln(£) - 3) - -i/; (^ ^^^3 ^ - tp ~ 


which, it is straightforward to show, is 0{£~‘^) as £ —>• 00 . Therefore the sum (^^ 1) convergent, and as it 

turns out, amenable to direct calculation. This is a somewhat involved and technical calculation, with details given 
in the Appendix. The result of this analysis is: 


1 00 ^ 

— Y^{2£+1) Pi = -l + 2zV31n 


e=o 


967r^ 


r(a) 

Lr(a)) 


- 12 —c ix,a) 
ax 


- {x,a) 


D-A 


D-3 




j=o 




iD-3) 


iD-3) 


+ {D - i){D - 5) {i; [{D -3){a+ 1)] +tP[l-{D- 3)a] - 21n(D - 3)} (. (2.64) 


This result has been confirmed numerically. 

By combining (I2.64|) with (12.59112.60p , and inserting the expression for a (j2.47p we finally arrive at a closed form 
expression for the renormalized vacuum polarization {(p^)ren of a massless conformally coupled scalar field, valid on 
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the horizon of the brane black hole 
D-3 


(0 )ren — 


487r2 


P-4)p-5) 




iP 


D-2 y'{D-A){D-2) 


2 V 3 


-iP 


D-2 y'{D-4){D-2) 

2 2^3 


-21n(L»-3) 


D-4 


167t‘^{D-3) 


X! -Cl+ 4) 


i =0 


1P 


6 j + 3p-2)- v^3(j^-4)(g-2) 
6(£i-3) 


+V' 


6j + 3(0 - 2) + v'3(i:'-4)(£)-2) 
6(£i-3) 


(2.65) 


We analyze this result, together with the renormalized 
vacuum polarization away from the horizon, in the next 
section. 


III. RESULTS 

In the previous section, we derived an exact, closed- 
form expression (12.651) for the renormalized vacuum po¬ 
larization on the event horizon of the brane black hole 
whose metric is given by (I 2 a . Outside the event hori¬ 
zon, the renormalized vacuum polarization is given as a 
sum of two terms, the first of which (I2.21all is given in 
closed form, while the second (I2.39P requires numerical 
computation. In this section we shall discuss our results 
for ((/)^)ren both On and outside the horizon. 

A. On the horizon 

We begin our discussion by examining the expression 
(12.651) for {(p^)ien on the event horizon. We see that for 
D = A and D = 5, the vacuum polarization on the hori¬ 
zon has the simple form 

which is in agreement with Candelas’ result for a 
Schwarzschild black hole space-time with D = A (4^ . 
Indeed, we may extend Candelas’ method to I? = 5 by 
noting that, for this case, (12.521) reduces to (2ttT = I 
(12.101) for D = 5 with rh = 11 [64|: 

Poe(r) = Pi(r), qoe(r) = Qe(r), (3.2) 

where Pi and Qi are Legendre functions. An application 
of Heine’s formula then allows us to calculate (12.651) , 
giving, for H = 5, 

('^^)ren= 2 ^’ (3.3) 


D 

(<^ )ren T — 1 

4 

l/(487r2) 

5 

l/(247r^) 

6 

0.0056676345 

7 

0.0065578987 

8 

0.0069605128 

9 

0.0069102842 

10 

0.0064254379 

11 

0.0055159981 


TABLE I. Values of the renormalized vacuum polarization 
(0">ren on the horizon ~ 1 oi a brane black hole for total 
number of space-time dimensions U = 4,..., 11. 


{<P\en 



FIG. 1. Renormalized vacuum polarization ((^^)ren on the 
horizon = 1 of a brane black hole for total number of 
space-time dimensions U = 4,..., 16. 


The on-horizon values of {(p‘^)ren for D = 4,..., 11 are 
given in Tab. U and the values for = 4,..., 16 are plot¬ 
ted in Fig. [H For D > 6 we give the numerical values of 
{(p'^)-cen in Tab. Uas they are more informative than the 
rather unwieldy analytic forms (I2.65|) . 

From inspecting Tab. |T] and Fig. [T] we see that {(p'^)ren 
reaches its maximum value on the horizon when the num¬ 
ber of bulk space-time dimensions D = 8, it then de¬ 
creases and first becomes negative at H = 15 (note it 
is still positive when D = 14). In fact we can see from 


which is in agreement with (EU- 
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direct consideration of (12.651) that the value of {(j)^)ren on 
the horizon continues to decrease toward —oo as D ^ oo. 
We may conclude from this that the magnitude of quan¬ 
tum effects near the horizon grows as the number of bulk 
space-time dimensions increases. We are working in units 
in which the event horizon radius is fixed to be equal to 
unity. In these units, the temperature (12.101) of the brane 
black hole increases linearly with the number of space- 
time dimensions D. Therefore, as H —>■ oo with the event 
horizon radius r/j fixed, the temperature of the black hole 
increases without bound and therefore the vacuum po¬ 
larization also increases in magnitude. One can therefore 
argue that as the number of extra bulk dimensions in¬ 
creases, the semi-classical approximation breaks down, at 
least in the vicinity of the black hole horizon. When we 
say that the semi-classical approximation breaks down, 
we mean that quantum effects are not small when D is 
large and the back-reaction of the quantum field on the 
space-time geometry can no longer be ignored. Eventu¬ 
ally a full theory of quantum gravity would be required to 
gain meaningful physical insights. In an upcoming paper 
[65l |. it will be shown that the same conclusion may be 
drawn for the vacuum polarization on the horizon in the 
bulk of the Schwarzschild-Tangherlini black-hole space- 
time (|2.1D . 


)analytic 

Next we turn to the renormalized vacuum polarization 
away from the horizon. The analytic part analytic 
(I2.21a|) . calculated using the metric function (12.2L takes 
the form 


(0 ) analytic — 


D-3 
1927r2/ _ 


Z)-3 


-2D D-5 


r,D-l 


,,2D-A 


(3.4) 


At the horizon r = I the expression ()3.4D simplifies to 


)analyti( 


r—1 


D-Z 
487r2 ■ 


(3.5) 


This matches the first term in the exact expression for the 
renormalized vacuum polarization on the horizon (12.651) . 
Furthermore, this is exactly {(jP')ren on the horizon when 
£) = 4 or 5 (13.11) . The remainder of (I2.65|) (which is 
nonzero only for D > 5) will come from the numeric part 
numeric- The expression (13.51) is positive for all D > A, 
and linearly increasing as the number of bulk dimensions 
increases. As r —oo, the expression (13.41) simplifies to 


) analytic 


(^-3)^ 

1927r2 


T2 
12 ’ 


(3.6) 


where T is the Hawking temperature (12.101) . The asymp¬ 
totic form (Irel) is simply the vacuum polarization for a 
quantum scalar field at temperature T in flat space-time. 
Like the renormalized vacuum polarization (13.51) at the 
horizon, the form ()3.6I) at infinity is also positive, but it 


)analytic 



FIG. 2. Analytic contribution (9 )anaiytic to the renormalized 
vacuum polarization on a brane black hole as a function of 
the radial co-ordinate r. The event horizon is located at r = 
1. The curves (from bottom to top) are for total number of 
space-time dimensions D = 4,..., 11. 


increases quadratically as the number of bulk dimensions 
increases. 

In Fig. [2] we plot ((^^) analytic (|3.4I) as a function of the 
radial co-ordinate r for total number of space-time di¬ 
mensions D = 4,..., 11 (from bottom to top curves in 
the figure). It can be seen that ((()^) analytic is positive ev¬ 
erywhere on and outside the event horizon, for all D. For 
H = 4, 5, 6 the maximum of (^^)anaiytic is on the horizon, 
and it is monotonically decreasing as r increases. For 
D >7 the maximum of analytic is a little outside the 
horizon. 

In [ 4 ^ it was found that analytic was the dominant 
contribution to the total renormalized vacuum polariza¬ 
tion {4>'^)ren On a foui'-dimensional Schwarzschild black 
hole. To see if this remains true on the brane when we 
have extra bulk dimensions, we now examine the numer¬ 
ical contribution (</>^)numeric- 


C. 




numeric 


The first step in calculating (())^)numeric is to find the 
mode functions Pui{r) and quii{r) by numerically inte¬ 
grating the radial equation (I2.14|) . The radial equation 
(12.141) has a regular singular point at the event horizon 
r = 1 and an irregular singular point as r —> 00. 

To find the mode functions Puje(r), which are regular 
at the horizon, we start our integration close to r = 1, 
using the power series 

OO 

Poje{r) = aj (r - 1)'"+^ , (3.7) 

j=o 

where v = ui/[D — Z) and Oq is set equal to unity, to give 
suitable initial values for Puie{r) and its derivative at the 
starting point. 
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The modes quji{r) are regular as r —>■ oo and have the 
following asymptotic series near infinity 

OO 

qu.t{r) ^ ^ (3.8) 

i=o 

where 

C = 

and we can set fop 
compute quitir) numerically, we found it easier to define 
a new function quii{r) = e^’^quiiij’) and numerically solve 
the resulting differential equation for qi^i{r), integrating 
from large r down to near the horizon. 

Since the mode sum in (j2.39|) involves cancellations be¬ 
tween the product of the mode functions Pcje{r), quji{r) 
and the WKB expansion terms Ckuitif), we require the 
mode functions to a high degree of precision. All numeri¬ 
cal calculations were performed in Mathematica, in which 
computations at the required precision are straightfor¬ 
ward to implement. The accuracy of our numerical in¬ 
tegration was checked by evaluating the normalization 
constant using (I2.15P at each value of r in our in¬ 
tegration grid. For all mode functions calculated, Cuj£ 
remained constant in r up to 26 significant figures. 

Once the mode functions have been found, the mode 
sum in (12.391) is straightforward to compute as the 
WKB terms (kuii (r ), while complicated, are algebraic ex¬ 
pressions. The numerical integrals Jk (12.381) are com¬ 
puted using standard numerical integration routines in 
Mathematica. 

Our results for ((^^) numeric on the brane as a function of 
the radial co-ordinate r are shown in Fig.|3]for total num¬ 
ber of space-time dimensions D = 4,11. In the top 
plot in Fig. [3] we show all eight curves for D = 4,..., 11 
(from top to bottom). In this plot, the curves for D = 4,5 
are virtually superimposed on the horizontal axis. To see 
these more clearly, in the bottom plot in Fig. [3] we show 
just the curves for = 4, 5. 

For £> = 4, in the bottom plot in Fig. [3] (9i^)numeric 
is positive everywhere outside the event horizon. We 
have compared our numerical values for ((^^) numeric with 
those tabulated in and the agreement is excellent. 
For D = 4, it can be shown that the numeric contri¬ 
bution (((>^) numeric vanishes on the event horizon r = 1. 
From Fig. it has a maximum just outside the event 
horizon and decays to zero quickly as r increases. Com¬ 
paring Figs. [3] and [21 for D = 4 the numeric contribution 
(</'^)numeric to the total renoriualized vacuum polarization 
(</'^)ren Is negligible compared to the analytic contribu¬ 
tion ((^^)anaiytic for all T, again in agreement with (d^ . 

For D >5, the numeric contribution (0^) numeric is neg¬ 
ative for all r > 1. When D = 5, comparing Figs. Inland 
[3] again reveals that (0^) numeric is negligible compared to 
(0^) analytic, but this is not the case lor D > 6. For all 


1 + w, 
1 , 


D = 4, 
D>5, 


(3.9) 


= 1 without loss of generality. To 




FIG. 3. Numeric contribution (0^)numeric to the renormal¬ 
ized vacuum polarization on a brane black hole as a function 
of the radial co-ordinate r. The event horizon is located at 
r = 1. (Top plot) The curves (from top to bottom) are for 
total number of space-time dimensions D = 4,..., 11. (Bot¬ 
tom plot) The curves (from top to bottom) are for D = 4, 5 
respectively. 


D > 5 we find that (0^)numeric has a minimum just out¬ 
side the event horizon, and increases toward zero as r 
increases. We also find that the magnitude of (0^)numeric 
for fixed r is increasing as the number of bulk space-time 
dimensions D increases. 


D. Total renormalized vacuum polarization 

We now combine the results of the previous two sub¬ 
sections to find the total renormalized vacuum polar¬ 
ization (0^)ren On the brane for a higher-dimensional 
Schwarzschild-Tangherlini black hole. Our results are 
presented in Fig.|4l where we plot (0^)ren as a function of 
the radial co-ordinate r for (from bottom to top curves) 
total number of space-time dimensions D = 4,..., 11. 
Our numerical calculations for (0^) numeric are valid only 
outside the horizon, for r > 1. We have extrapolated our 
numerical results for (0^)ren for r > 1 to the horizon at 
r = 1 and find excellent agreement with the exact results 
on the horizon given in Tab. [H 

For = 4,..., 11 we see from Fig.|4|that the renormal¬ 
ized vacuum polarization (0^)ren is positive everywhere 
on and outside the event horizon. For H = 4, 5, 6 it has 
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<0'>ren 



FIG. 4. Total renormalized vacuum polarization ((^^)ren on 
a brane black hole as a function of the radial co-ordinate r. 
The event horizon is located at r = 1. The curves (from 
bottom to top) are for total number of space-time dimensions 

its maximum value on the horizon, and is monotonically 
decreasing as r increases. For D > 7 we find that ((^^)ren 
has a minimum value on the horizon and is monotoni¬ 
cally increasing as r increases and we move away from 
the event horizon. For all the numeric contribution 
(</'^)numeric becomes insignificant for large r, and the ana¬ 
lytic contribution analytic dominates. However, close 
to the horizon the numeric part ((^^) numeric makes a sig¬ 
nificant contribution to the total. 

Far from the event horizon of the brane black hole, 
the vacuum polarization ((^^)ren increases rapidly as the 
number of bulk space-time dimensions D increases. As 
explained earlier, this is to be expected because we are 
working in units in which the event horizon radius is 
fixed to be unity, and in these units the temperature of 
the brane black hole (12.101) increases linearly with the 
number of bulk space-time dimensions D. 

Close to the horizon, due to the significant nega¬ 
tive contribution to the total from the numeric part 
numeric, the vacuum polarization is rather smaller 
than it is far from the black hole. As shown in Fig.[T] for 
H = 4,..., 11 the vacuum polarization is always positive 
on the horizon and has a maximum when D = 8. 


IV. CONCLUSIONS 

In this paper we have calculated the renormalized vac¬ 
uum polarization ((/)^)ren for a massless, conformally cou¬ 
pled, quantum scalar field 0 propagating on the brane 
of a higher-dimensional Schwarzschild-Tangherlini black 
hole. Our four-dimensional, on-brane metric, is static 
and spherically symmetric and corresponds to a slice of 
the higher dimensional black hole geometry. The metric 
functions depend on the total number of space-time di¬ 
mensions D, and in this paper we focus on I? = 4,..., 11, 
in common with the literature on Hawking radiation from 


brane black holes [1, [l^, flbl - l^ . The Hawking radiation 
of a scalar field from such a brane black hole has been 
studied in depth (see, for example, [s^) but the Hawking 
fluxes can be computed without recourse to renormaliza¬ 
tion. To the best of our knowledge, we have presented in 
this paper the first computation of an expectation value 
requiring renormalization on a brane black hole space- 
time. 

Since the on-brane metric is four-dimensional, static 
and spherically symmetric, we have been able to employ 
well-established methodology [l^, [H, |4l|, to com¬ 
pute the renormalized vacuum polarization. We work 
on Euclidean space-time and consider the scalar field to 
be in the Hartle-Hawking state [s^, a thermal quan¬ 
tum state at the Hawking temperature T. We use co¬ 
variant geodesic point separation to regularize the Eu¬ 
clidean Green’s function and then the divergent subtrac¬ 
tion terms are given by the Hadamard form. Away from 
the horizon, the space-time points are split in the tem¬ 
poral direction. In this case the renormalized vacuum 
polarization can be written as the sum of two terms, the 
first of which is a simple closed-form expression, while the 
second requires numerical computation. On the horizon, 
radial point-splitting is used and an exact (but compli¬ 
cated) closed-form expression for the renormalized vac¬ 
uum polarization is found. 

For D = 4,..., 11 we find that the renormalized vac¬ 
uum polarization is positive everywhere on and outside 
the event horizon. Far away from the horizon, its value 
increases rapidly as the number of bulk space-time di¬ 
mensions D increases. This is because the on-brane 
temperature of the black hole increases linearly with in¬ 
creasing D for fixed event horizon radius Vh and we are 
considering a thermal quantum state. Ultimately, for D 
very large (with fixed r^) the semi-classical approxima¬ 
tion used here breaks down in the sense that quantum 
effects are no longer small and the back-reaction of the 
quantum field on the black hole geometry can no longer 
be ignored to first order. This is also the case for very 
small black holes which have entered the Planck phase 
of their evolution Q. In either of these regimes {D very 
large or very small black holes) a full theory of quantum 
gravity is required to model the behaviour of the black 
holes. Close to the horizon, the renormalized vacuum po¬ 
larization becomes negative when U > 15 and decreases 
rapidly as D increases further. This again ultimately 
leads to a breakdown in the semi-classical approximation 
when D is very large. 

In this paper we have taken the quantum scalar field 
to be in the Hartle-Hawking state since this state is 
the easiest in which to compute renormalized expectation 
values. The quantum state of interest for simulations of 
brane black holes at the LHC [3, [T^ is the Unruh state 
since this state represents an evaporating black hole. 
Since differences in expectation values between two quan¬ 
tum states do not require renormalization, it would be 
comparatively straightforward to compute the renormal¬ 
ized vacuum polarization in the Unruh state using the 
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results presented here for the Hartle-Hawking state. 

We have considered the renormalized vacuum polar¬ 
ization ((^^)ren becausc it is the simplest nontrivial ex¬ 
pectation value for a quantum scalar held. While it is 
a scalar (and hence cannot, for example, distinguish be¬ 
tween the future and past event horizons of a black hole), 
it nonetheless shares some physical features with the 
renormalized stress-energy tensor (T^[/)ren. The renor¬ 
malized stress-energy tensor is of particular interest be¬ 
cause it governs the back-reaction of the quantum held on 
the space-time geometry via the semi-classical Einstein 
equations (HU). However, since (T)ji/)ren is a tensor ob¬ 
ject and involves derivatives of the Green’s function, it is 
more complicated to calculate than {(fP),:en, although the 
standard methodology [H, should be applicable for 
the on-brane metric considered in this paper. It would 
be interesting to investigate whether the effects we have 
found here for the renormalized vacuum polarization (in¬ 
creasing magnitude with increasing bulk dimension D 
and negative values near the horizon for large D) are 
present also in (T)ji/)ren- 

There are other extensions to this work on the brane 
which would be of interest. In this paper we have con¬ 
sidered a massless, conformally coupled scalar field. It 
is known that including a mass suppresses the Hawk¬ 
ing radiation on the brane and so the emission of light 
quantum fields is of the greatest interest phenomenolog¬ 
ically. For a four-dimensional Schwarzschild black hole, 
the Hawking radiation is dominated by scalar fields but 
as the number of bulk dimensions increases, the on-brane 
emission of particles with non-zero spin becomes of com¬ 
parable magnitude to scalar emission [s^ . Emission of 
spin-half quanta has particular phenomenological impor¬ 
tance given the large number of fermion degrees of free¬ 
dom in the Standard Model. Therefore the computation 
of (T^i/)ren for non-zero spin fields, while technically chal¬ 
lenging, would be of interest. 

Here we have restricted our attention to static, nonro¬ 
tating black holes. Hawking radiation on the brane from 
rotating black holes has been extensively studied (see, 
for example, [l^, [2^ for reviews) but the computation 
of renormalized expectation values everywhere outside 
the horizon of a four-dimensional rotating black hole has 
proved intractable to date. 

Finally we emphasize that our results in this paper are 
for the renormalized vacuum polarization on the four¬ 
dimensional brane metric of a higher-dimensional black 
hole. For smaller values of D, emission on the brane dom¬ 
inates that in the bulk [^, , but as the number of 

bulk dimensions D increases, Hawking radiation of gravi¬ 
tons (and possibly scalar fields) in the bulk does become 
significant. While the formalism underlying the renor¬ 
malization of expectation values for quantum scalar fields 
on higher-dimensional space-times is known [3^ . [^ . a 
methodology for detailed computations in the full space- 
time exterior to the event horizon has yet to be devel¬ 
oped. We leave these questions for future work. 


Appendix A: fie summation 

In this Appendix we give the details of the derivation 
of dam). We begin with 

(3e= i In (^ -I- a) -f i In (£ -I- a) - ln(T) - 3) 




D-3 


D-3 


-I- a ) — -'ip 


i + l 
D-3 


(Al) 


In order to find the sum (2^ + 1) we consider 

instead the sum 






(A2) 


e=o 


Our strategy is to expand E(z) and its derivative about 
2 = 1, which we may then combine to give 


5^(2£+1)/3^ = 2S'(1) + S(1). 


(A3) 


e=o 


upon taking the limit z —>■ 1. 

We begin by examining the expression 


e=o 




+ q z 


-E 

i—p 


ip 


q]z^-'tp(-+q-l]z 


.(A4) 


Introducing a new variable k = i — p and making use of 
Eq. (8.365.1) of 


tpiv -|- 1) — 'ip{v) = 


1 


(A5) 


we arrive at 


e=o 


E^(-+9)^^ = 


•p-1 


1- zP 


to VP ' 


k +pq 


(A 6 ) 


From Eq. (8.365.6) in [^, we have 

^■ip i- + zj =p [ipipz) - ln(p)], (A7) 

e=o 2 

and so the right hand side of (IA6I) can be written as 

1 


e=o 


1- zP 

+p{z’’j2 




OO 1 . 

ZK 


e=o 


k + pq 


■ '0 ipq) - ln(p) 


■ (AS) 
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The definition of the hypergeometric function gives, 

OO 1. 

1. K ^ Z 

— 2 F 1 {l,pq-,l+pq]z) = > 

Tin 


pq 


k=0 


k +pq 


(A9) 


into three separate parts 

E(z)=^1+52+^3, 


(AlO) 


Combining the above work, we may split our sum 'S{z) where 

I 


S'! = i ^ [ln(£ + a) + ln(t' + a)] 


e=Q 

' ' i=o 


5 , = - 


D-3 


D-3 

D-3 


+ a + 1 ] + 


e + i 

D-3 


— a 




2(l-z-D-3) \i^D-3){a + l) 
■.d-3 


[l-[D- 3)a] + 


l-[D- 3)( 


2F1 (1, {D - 3)(a + 1); 1 + (i^ - 3)(a + 1); z) + V' [{D - 3){a + 1)] 

,M - (!> - 3)0; 2-in- 3)0; .)l - ffcl) - =« 


1-z 


1 - z 


D-3 


(All) 


We now proceed to expand about z = 1, retaining the 0(z — 1) terms as we will eventually be taking the derivative 
of the resulting expressions with respect to z. We will deal with S 2 and S 3 first as they are the most straightforward 
to work with. The required expansions are readily obtained using Mathematica; 


^2 = 


1 


2(D - 3) 


D-4 


i=0 


if) 


D-3 


+ a + 1 + tp 


z- 1 
'4{D-3) 


' D-A 


Y^I{^ + 3-D) 


t=o 


V' 


D-3 


l+l 

D-3 


+ a + 1 + '0 


D-3 


+ 0[z-lf 


(A12) 


and 


83 = - {{D - 3)aiP [{D - 3)(a + 1)] - [{D - 3)a + £> - 4] 0 [1 - (£> - 3)a] + (£> - 4) [ln(D - 3) - 1]} 
+ ^ {[{D - ^)iD - 5) - 6iD - 3)a] 0 [{D - 3)(a + 1)] 

+ [(£> - 4)(£i + 1) + 6{D - 3)a] 0 [1 - (£> - 3)a] + 6{D - 4) - 2{D - 2){D - 4) ln(£i - 3)} 

- i°00;+'' +o[(.-if]. 

Returning to ^i, from Eq. (9.550) in we have that 

.1 


E 

t=3 


(v + £)^ 


= <^(z,s,v), 


1^1 < 1 , 


where <l’(z,s,r;) is Lerch’s transcendent function. Therefore 


^ z^ In (£ + u) = s, d)|s=o 




(A13) 


(AM) 


(A15) 


and so we can express Si in the form 




ds 


$(z,s,a) 


s=0 


as 


s=0J 


(AI6) 


To expand this expression about z = 1 we make use of the following result for $(z, s,z;) (see Sec. 1.11 in |68l |') 

n s -1 00 


4>(z,s,u) = 


r(i-s) 


In 


+ z - k,v) 


k—0 


{InzY 


(A17) 
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Taking derivatives of both sides with respect to s and setting s = 0 yields 


A. 

ds 


5>(z,s,t;)U=o = ^ |ln(r(r;)) -ln(v^) + ^ k, v) 


s=0 


(Inz)'^ 
kl 


where we have used the result 


(A18) 


ds 


C{s,v) 


= ln(r(u)) — ln(\/^). 


s=0 


Expanding (IA18I) about z = 1 and inserting the resulting expression into (IA16I) yields: 

S'! = — [ln(27r) — 1 — Inr(a) — Inr(a)] H-~ 61n(27r) + 12 [alnr(a) + alnr(a)] 

+ i"(i;Ji+A o[(.-if]- 




s=0 


+ ^C(»-i,a) 


s=0 


Combining (jA12[ lAT^ IA20I) yields 


E(z)=Fi+F2(z-1) + 0[(z-1)2] 


where 


1 


In 


27r 

D-3 


+ {D-3) [\n{D - 3) - 1] - Inr(a) - Inr(a) 


+ aiD-3)'iP[iD-3)ia + l)] 


- [{D - 3)a +{D- 4)] 3)a] } + 


1 


D-i 


2{D-3) 






D-3 


+ a + 1 


■'0 


£+1 


D-3 


1 


F 2 = —\6D- 19 + 21n(H- 3) -61n(27r) - 2(D - 3)^ ln(£» - 3) + 12alnr(a) + 12alnr(a) 


+ [{D - 5)iD - 4) - 6{D - 3)a] 0 [{D - 3)(a + 1)] 

+ [{D - 4)(H - 5) + 6{D - 3)a] 0 [1 - (D - 3)a] - 12 

D-4 


ds 


C(s - l,a) 


s =0 


+ ^C(»-i,a) 


s=0. 


+I775^E^(« + 3-i>) 


4(D - 3) 


e=o 


0 


D-3 


+ a + 1 + 0 


D-3 


— a 


(A19) 


(A20) 


(A21) 


(A22) 


Finally, by taking the derivative of (|A21ll with respect to 
2 and setting ^ = 1 we obtain: 


^(2£+1)/3£ = 2F2 + Ei, 


(A23) 


i=0 


which simplifies to the expression contained in (12.641) . 
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